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Abstract 

Continuum mechanics (e.g., hydrodynamics, elasticity theory) is based on the assumption that a 
smah set of fields provides a closed description on large space and time scales. Conditions governing 
the choice for these fields are discussed in the context of granular fluids and multi-component fluids. 
In the first case, the relevance of temperature or energy as a hydrodynamic field is justified. For 
mixtures, the use of a total temperature and single flow velocity is compared with the use of 
multiple species temperatures and velocities. 

PACS numbers: 45.70.-n,47.35.-i,45.70.Nd,51.10.+y 
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I. INTRODUCTION 



Continuum mechanics (e.g., hydrodynamics, elasticity theory) provides the macroscopic 
description for a wide class of physical systems and states in terms of a few space-time 
fields. The derivation and justification of the equations for these fields from the underlying 
Newtonian or quantum mechanics is a well- studied and open problem of non-equilibrium 
statistical mechanics [l-3|. The central conceptual and practical problem is the reduction of 
the many degrees of freedom in the microscopic description to a closed description in terms 
of a few chosen macroscopic fields. Qualitatively, this is understood in terms of restrictions 
on the conditions for applicability of continuum mechanics. In particular, it is expected 
that such a reduced description can apply only for slowly varying fields that dominate the 
description on sufficiently long time scales. Under such conditions, all other fast degrees of 
freedom have decayed to zero leaving the possibility for a reduced description. 

The choice of fields is therefore an important first step in any derivation of continuum 
mechanics. This in turn requires understanding relevant time scales active in the system. 
For a simple molecular fiuid, there is a time independent uniform state, equilibrium, which 
is the reference state for the dynamics. The longest time scale for approach to equilibrium 
is set by the relative size of the spatial deviations of the relevant macroscopic fields y, 

~ Ay /lie, where Ay ~ ?/(r + ro,t) — y{r,t), y^ is the equilibrium value of y, and tq 
is a distance of the order of the mean free path. In the following, /i is referred to as 
the uniformity parameter. The local conserved fields of number, energy, and momentum 
densities obey balance equations expressing their time derivatives as gradients of fiuxes. 
These derivatives are therefore proportional to the uniformity parameter, identifying a time 
scale that diverges as the system approaches equilibrium. A necessary condition for choosing 
a set of fields, therefore, is that they should include all local conserved densities. In the case 
of broken symmetries (e.g. solids, liquid crystals, superfiuids) there are additional non- 
conserved fields whose time derivatives are proportional to the uniformity parameter, and 
therefore must also be included in the set of chosen fields j^, 5|. Finally, there are systems 
with internal processes on long time scales (polymers, plasmas) that remain finite in the 
uniform limit. Their associated fields can be included in the chosen set if these time scales 
are long compared to all other non-conserved excitations. 

In the following, attention will be restricted to the fiuid phase, for which the macroscopic 
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description is hydrodynamics. Also it is assumed that there are no slow excitations associ- 
ated with a broken symmetry. For reasons just noted, the hydrodynamic fields must include 
all conserved densities as a necessary condition to be complete on the longest time scale. The 
question addressed here is whether this minimal set can be extended to include additional 
fields, and if so, whether it is appropriate to do it. This issue has been discussed extensively 
in the context of extended irreversible thermodynamics (EIT) j^. In EIT, additional fields, 
such as the fiuxes in the balance equations for the conserved densities are included as inde- 
pendent hydrodynamic fields. One motivation there is to resolve paradoxes of Navier-Stokes 
hydrodynamics, usually associated with a misuse of those equations for conditions where 
they do not apply (e.g., initial slip, boundary layers, shock fronts). Another approach using 
additional fields is the moment method of Grad [7]. Here, motivated by current applications 
of hydrodynamics to granular fiuids [s], related but different issues are addressed. 

Section [TTl outlines the basic ingredients for a derivation of hydrodynamic equations near 
a state of uniformity, from the underlying equations of nonequilibrium statistical mechanics 
9|. It relies on exact balance equations for the chosen fields, obtained from the Liouville 
equation, and on a special "normal" solution to that equation. Construction of a approx- 
imate normal solution is obtained by expansion in the uniformity parameter, leading to 
Navier-Stokes level hydrodynamics. The derivation does not make any explicit limitation on 
the chosen fields, other than the existence of a uniform reference state expressible in terms 
of those fields. However, it does not provide the context under which such a solution is 
applicable. It is argued that a necessary condition must be sufficiently long time scales. 

In Sec. mil the example of a simple granular fiuid is considered. The uniform state 
in this case is not equilibrium, but one for which the temperature decreases monotoni- 
cally due to inelastic collisions. The most significant differences from a molecular fluid are 
non-conservation of energy and a uniform state that is time dependent. The number and 
momentum densities are still conserved, so they are included in the set of hydrodynamic 
flelds. The question of including the non-conserved energy density (or, equivalently, the 
granular temperature) is discussed. It is concluded that the minimal set of hydrodynamic 
flelds must include the temperature. 

In Sec. IIVI binary mixtures are addressed, for both molecular and granular fluids. Two 
sets of flelds are considered for hydrodynamic descriptions. One is the conserved number 
densities, flow velocity associated with the conserved total momentum, and the temperature 
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associated with the total energy. The second one is an expanded set consisting of the 
conserved number densities, the species flow velocities associated with the non-conserved 
species momenta, and the species temperatures. It is argued that the second description, 
while more detailed, has no predictive value beyond the simpler hydrodynamic description 
on the relevant large space and time scales. The presentation is summarized in the last 
section with some concluding remarks. 

II. NORMAL SOLUTIONS, CONSTITUTIVE EQUATIONS, AND HYDRODY- 
NAMICS 

The most general notion of a hydrodynamic description is a closed set of equations for 
a set of hydrodynamic fields, denoted by {i/a}- The terminology "closed" means that they 
obey a set of equations of the form 

dtya{r,t) = N^{r,t\{yp}), (1) 

where the space, r, and time, t, dependence of the right side is entirely determined from 
the fields themselves. In general, A'"Q,(r, t| {y^}) is a nonlinear functional of the fields. The 
derivation of such equations procedes in two steps. First, balance equations for the fields 
are obtained directly from the Liouville equation, 

%a(r,t) + V-j,(r,t) = s„(r,t). (2) 

Here, the fluxes jo, (r, t) and sources Sa{r,t) are linear, time-independent functionals of the 
solution to the Liouville equation. In general, these functionals cannot be characterized en- 
tirely by the hydrodynamic fields, so although ([2]) are exact equations they do not constitute 
a closed set of equations for the fields in terms of themselves. This closure occurs for the 
special conditions of a "normal" state, defined as a solution to the Liouville equation of the 
form 

Pn (r, t) = Pn (r| = pn {{qij} , {v^} , (q^, t)}) , (3) 

where Qi and Vi are the position and velocity of particle i, and Qij = qi — qj. For such solu- 
tions, all time dependence and the breaking of translational invariance occur only through 
the hydrodynamic fields. The notation in the second equality indicates that the functional 
of the fields pn (F | {tja}) is equivalent to a function of the fields at the particle positions of 
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the specified T point. For a normal solution, the definitions of {ja {r,t)} and {sa{r,t)} as 
linear functionals of this solution lead to "constitutive equations" 

ia{r,t) = ia{r,t\{yp}), s„(r, t) = s«(r,t| {y^}). (4) 

With the fiuxes and sources now determined by the fields, the balance equations (j2]) become 
closed hydrodynamic equations of the form ([1]). 

Since the balance equations are exact consequences of the microscopic dynamics, ques- 
tions about the existence of a hydrodynamic description turn on the necessary conditions 
for a normal solution (e.g., appropriate choice for the set of fields {ya}, restrictions on the 
space and time scale). Substitution of ([3]) into the Liouville equation gives the equation a 
normal solution must satisfy for a given choice of fields {ya}, 

I . t\ {yp}) + Lp„ = 0, (5) 

where summation over repeated indices is implicit, L is the Liouville operator of the system, 
and Na{r,t\ {y/3}) is the nonlinear functional in the hydrodynamic equations ([T]), defined 
in terms of the fiuxes and sources. The latter in turn must be determined self-consistently 
in the form of constitutive equations (jl]), from the solution to The resulting solution 
is implicit, being a function of the phase point F and a functional of the unknown fields. 
The space and time dependence of these fields is then obtained from a solution to the 
hydrodynamic equations for given initial and boundary conditions, completing specification 
of the normal solution. 

Determination of a normal solution in this way is quite difficult, but its approximate 
construction can be carried out explicitly for states with small uniformity parameter n 
(described above), by expansion in this parameter. This construction of the normal solution 
by expansion around the corresponding homogeneous solution will be referred to as the 



Chapman- Enskog method [10|, This presumes the existence and determination of a 
reference homogeneous solution ph to Eq. (I5l), 

^^ySa{{yf^m) + Lp, = 0, (6) 

together with the homogeneous hydrodynamic equations 

dtya{t) = sa{{yp{t)}). (7) 
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The familiar example is a one component atomic fluid near spatially homogeneous equilib- 
rium. If the flelds chosen are all associated with local conserved densities, then the sources 
Sa vanish and (jS]) becomes the condition for a stationary equilibrium state. More generally, 
if the flelds have sources, the homogeneous solution depends on the time dependence of these 
flelds generated by the sources. Effects of small spatial gradients are described by returning 
to ([5]) and expanding the solution about ph to linear order in p, 

Pn (r I {y^}) = Ph (F; {y^}) + A„(r; {y^ (r, t)}) ■ pVy^ir, t) + . . . (8) 

The dots indicate terms of higher order in p, and ([5]) provides a linear inhomogeneous equa- 
tion for determination of Aq,. With this result, the constitutive equations also are determined 
to flrst order in the gradients, providing Na and the hydrodynamic equations to this order. 

If the flelds are chosen to be the local conserved densities of mass, energy, and momentum, 
then the sources Sa{r, t) all vanish, and the fluxes ja(r, t) are proportional to the uniformity 
parameter. Construction of the normal solution by the Chapman-Enskog method determines 
the time derivatives to be proportional to the uniformity parameter as well. The leading 
approximation yields the perfect fluid Euler equations, and the next order gives the usual 
Navier-Stokes equations. The derivation provides their context for applicability as well in 
this case: states must be near uniform equilibrium, with small fractional variations in the 
flelds over a mean free path, and space and time scales large compared to the mean free 
path and mean free time, respectively. 

It would appear that this approach is applicable for any choice of flelds, leading to many 
forms for the macroscopic description of a system. However, the conditions under which 
a normal solution could be expected must be considered further. Consider flrst a normal 
fluid with elastic collisions in an initial non-equilibrium state with specifled hydrodynamic 
flelds {ya (r, t = 0)}, whose values vary smoothly across the system. In each small region 
of dimension larger than a mean free path, the phase space density p{T,t) approaches a 
Gibbs density characterized by the hydrodynamic flelds at its central point r. Subsequently, 
exchange of mass, energy, and momentum tends to equilibrate these flelds to uniform values 
(or to steady values if the system is driven). The flrst stage, approach to a universal 
form for the velocity distribution, occurs after a few collisions. This establishes the normal 
form of the solution, where the hydrodynamic flelds and their gradients characterize the 
state. Deviations from the equilibrium distribution are due to fluxes of mass, momentum. 



6 



and energy across the cells. These fluxes are proportional to the differences in values of 
the fields (i.e., to their spatial gradients). The second stage is the slower evolution of the 
distribution through the changing values of the fields, according to the "usual" Navier-Stokes 
hydrodynamic equations. In this scenario, there is a characteristic microscopic time scale 
(mean free time) and a characteristic microscopic space scale (mean free path) that set the 
limits for validity of the normal solution. Beyond these limits, the dynamics is complex and 
not captured entirely by that of the hydrodynamic fields. 

Extending the hydrodynamic fields to include non-conserved fields necessarily introduces 
additional time scales in the hydrodynamic solutions, due to the sources s„ associated with 
the new fields. Thus the hydrodynamics of the second stage now has variations on the 
same time scale as the first stage. Furthermore, the additional fields give only a partial 
description of the short time scale dynamics, and therefore do not extend the domain of 
applicability for the hydrodynamic description. The qualitative initial stage described above 
for a general initial preparation of the system is still required for any reduction in the large 
number of degrees of freedom active in the system to decay. The hydrodynamics can be 
closed (Markovian) only after such a reduction. Even during the second stage, only the 
slower dynamical processess can be described by hydrodynamics, as there are continual 
fluctuations of the fast degrees of freedom that would prohibit isolation of the corresponding 
fast hydrodynamic modes. In summary, the choice of fields for hydrodynamic equations is 
somewhat arbitrary, but their dominance on some time scale is a central property for utility 
as a macroscopic description of the system. Predictions on shorter time scales are difficult 
to correlate with measurements, and potentially obscuring the simplicity of the larger time 
scale dynamics. On the common long time scale required for validity, it would appear that 
the normal solutions should be equivalent, with the more detailed fields becoming functionals 
of the simpler fields and being "slaved" to them. 

Exceptions to this conclusion occur when the uniform reference state is not the equilib- 
rium state, and may have an inherent dynamics of its own. This is the case for granular fiuids. 
In the following, the hydrodynamic descriptions for two examples using non-conserved fields 
are discussed: a simple one-component granular fluid for which there is an energy source, 
and binary mixtures where non-conserved flelds of individual species may be of interest. 
Attention is limited to states near a corresponding universal homogeneous state (homoge- 
neous cooling or equilibrium state) for simplicity. The question addressed is "what are the 
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appropriate Navier-Stokes hydro dynamic equations for these more complex conditions, and 
under what conditions do they apply?" . 



III. GRANULAR FLUID 

A one-component granular fluid consists of "particles" (grains) of mass m with short 
ranged collisions that conserve number and momentum, but not energy. By analogy with 
a molecular fluid it is reasonable to expect a hydrodynamic description in terms of the 
local number, momentum, and energy densities. It is usual to replace the momentum by 
its associated flow velocity u, and the energy by a granular temperature T. This is just a 



change of variables. The exact macroscopic balance equations for n, u, and T are 12 1 



DtTi + nV • u = 0, (9) 
Au+^V-P = 0, (10) 

AT + 1^ (V ■ q + P : Vu) = -CT. (11) 

where Df = dt + u ■ V is the material derivative, and p = nm. For simplicity, the particles 
have been represented as hard spheres for which the temperature is defined in terms of the 
energy density e as e = 3nT/2 + pu^/2. The heat flux q and pressure tensor P are defined 
as the energy and momentum fluxes in the local rest frame of a fluid element moving with 
velocity u. They are related with the energy and momentum fluxes in the laboratory frame, 
s and T, by 

q = s— eu — u ■ P, P = T — puu. (12) 

These balance equations have the same form as for a molecular fluid, except for the presence 
of a source —(T in the temperature equation, resulting from the inelasticity of collisions 
among the particles. The parameter ( will be referred to as the cooling rate. Application of 
the Chapman-Enskog method for a normal solution to first order in the uniformity parameter 
leads to constitutive equations for the heat and momentum fluxes, q and P, similar to 
Fourier's law and Newton's viscosity law jl3| 

q = - AVT - fiVn, Pij = p6ij - t] (ViUj + VjUi - ^^ijV • - uSijV ■ u. (13) 
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In addition, the constitutive equation for the coohng rate ( up to first order in the gradients 
is given by 

C = Co + CiV ■ u. (14) 

The pressure p and the transport coefficients in these expressions, A, /i, rj, k, (q, and (i, are 
also determined from the normal solution as functions of n and T. 

Equations fl9|)-f lT4|) constitute the Navier-Stokes order hydrodynamic equations for a sim- 
ple granular ffuid. In contrast to a molecular ffuid, solutions to the granular hydrodynamic 
equations have two time scales, one set by the uniformity parameter and another set by the 
energy loss rate r ~ (q^. The former becomes large as the system approaches uniformity, 
while the latter remains finite in uniform systems. From the discussion of the last section, 
in order for the hydrodynamic description to dominate the complex microdynamics, both of 
these time scales should be large compared to the mean free time. In practice this appears 
to be the case, at least if the degree of inelasticity is not too large (since Co vanishes in the 
elastic limit). 

As the system approaches uniformity, the two hydrodynamic time scales become very 
different and the question arises of whether the temperature is still a relevant hydrodynamic 
field. It is a "fast variable" on the time scale for the spatial variations, so perhaps there is a 
simpler description in terms of n and u alone. To see that this is not the case consider Eq. 
(|6]) for the reference homogeneous normal solution, 

(-CT^ + p. = 0. (15) 

To exclude the temperature as a relevant variable, ph should be independent of the temper- 
ature, i.e. T should represent a transient dynamics that vanishes or becomes constant on a 
sufficiently long time scale. Then ( lT5l) would become Lph = for the reference homogeneous 
normal solution on this time scale. However, there is no solution to this equation for finite 
inelasticity. This can be easily seen by considering the average rate of change in the energy 
E, 

dt{E) = - j dTELpj,. (16) 

This can never vanish since the coUisional cooling continues as the system evolves. Thus 
the temperature is an inherent property of all solutions and cannot be neglected on any 
time scale. Although (o characterizes a rapid variation relative to the spatial time scale. 
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it persists as a modulation of that latter slow scale. The reference uniform state has a 
dynamics (cooling) that is inherited by spatial deviations from that state and the system 
continues to cool rapidly even after it has "aged" an arbitrarily long time. 

It might be argued that for special boundary conditions the hydrodynamic equations 
would support a solution with stationary temperature, such that variations of the tem- 
perature are indeed transient and rapidly approaching its stationary value. Subsequently, 
the residual relaxation of spatial variation might be described by n and u alone, and only 
parameterized by a constant temperature. While such a contrived description might be 
possible on a case-by-case basis, this does not constitute an alternative hydrodynamics. A 
hydrodynamic description is a set of closed equations for chosen fields that describes a class 
of states, e.g. Navicr-Stokcs equations for states near homogeneity, whose form is universal 
and independent of specific initial or boundary conditions. 

The granular fluid with the source —(^T can be contrasted with a molecular fluid including 
non-conserved flelds. In the latter case, the sources vanish for the homogeneous (equilibrium) 
state so that solutions to Lpe = do exist. There is then no restriction on the removal of 
the non-conserved fields in that case, for a simpler hydrodynamics. In contrast, the granular 
source discussed here docs not vanish for its homogeneous state and hence the dynamics due 
to this source must be included in any set of hydrodynamic fields. 

In summary, the Navier-Stokes hydrodynamics for a simple granular fluid is given in 
terms of n, u, and T just as for a molecular fluid. Due to the inelasticity of coUisions, there 
is a new time scale set by the cooling rate that does not increase as the system approaches 
spatial homogeneity. As for molecular fluids, the normal solution and hydrodynamics applies 
for times long compared to the mean free time when the transient dynamics of microscopic 
degrees of freedom has decayed and only that of the fields remains. However, the granular 
hydrodynamics retains a modulation of the spatial relaxation due to continual cooling of the 
system. 

IV. BINARY MIXTURES 

Consider next a granular or molecular fluid composed of two different species (e.g., dif- 
ferent mass, size or collision law). If there are no reactions or fragmentation, the number of 
particles of each species is conserved and their local densities, rii, are appropriate indepen- 
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dent hydrodynamic fields. Other properties of interest might be the species temperatures, 
Tj, and the species flow velocities, Uj. However, these are not associated with any conserved 
quantity even for a molecular fluid (there is continual energy and momentum exchange be- 
tween species, even though the total energy and momentum are conserved). In this section, 
the possibility of including species properties, beyond their number densities, as hydrody- 
namic fields is discussed. 

The species numbers and flow fields, related to the corresponding total fields 

n and u by 

ni + ^2, pu =nimiUi-Fn2m2U2, (17) 

where the mass density is p = mini + 7712^2 — Pi + p2- The temperature is determined 
from the total energies which are decomposed into their kinetic temperature and convective 
energies according to 

3 1 3 1 

e = -nT -pu^ = ei 62, = -UiTi + -piuf. (18) 

Two possible hydrodynamic descriptions are considered here, one for the six fields {rij, u, T}, 
and another for the ten fields {n^, Uj, Tj}. 



A. Single temperature, single velocity balance equations 

The exact macroscopic balance equations for ni, n2, u, and T are 

Dtni + mW -u + W -ii^O, (19) 
Au + ^V-P = 0, (20) 
AT--V- ^j, + |-(V-q+P: Vu) = -CT, (21) 

^ 1=1,2 ^ 

where ji is the number flow for species i relative to the local flow, jj = rij (uj — u). To 
obtain hydrodynamic equations for this choice of fields, constitutive equations must be 
obtained from the corresponding normal solution to give jj, q, P, and C as functionals of 
these fields. As with the one component granular fluid, there are two time scales associated 
with this choice of hydrodynamic fields, that associated with the uniformity parameter and 
that associated with the cooling rate. Both should be larger than the mean free time. 
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B. Two temperatures, two velocities balance equations 



The balance equations for the extended set of ten fields have the form 







(22) 



PiDitUi + V ■ Pi 
(V ■ qi + Pi : VUi) 



A,. 



(23) 



DuTi + — 




Ui ■ A,; 



■J? 



(24) 



where now Da = 9j + Uj ■ V, and the heat flux and pressure tensor Pj are now defined 
as the energy and momentum fluxes in the local rest frame of a fluid element moving with 
velocity Uj (compare with Eq. f|T2|) ). The sources Aj and Q depend on rij, Tj, and Uj. The 
normal solution for this choice of fields is needed to get the constitutive equations for q^, Pj, 
Aj, and Q as functionals of these fields. 

Note that q 7^ qi + q2 and P 7^ Pi + P2, in general. The balance equations for ni,n2, u, 
and T, Eqs. ( 1T9|) -(I2T]). follow from Eqs. fl22l)-(l24l) only if the sources have the properties 



Here Sj and Tj are the energy and momentum fluxes of species i in the laboratory frame. 
These relations are proved in Appendix |X1 where explicit expressions for the sources are 
obtained. 

The coefficients Aj are proportional to the collision frequency, and hence introduce dy- 
namics on the short mean free time scale. The rates Q have two contributions. One is 
positive and proportional to the degree of dissipation, like ( in Eq. ( 12T|) . However, there is 
an additional contribution of indefinite sign and which does not vanish in the elastic limit. 
It represents the collisional transfer of energy between different species. This latter contri- 
bution is of the order of the collision rate and consequently generates dynamics on the mean 
free time scale. There are now three hydrodynamic time scales: the largest one set by the 
uniformity parameter, that set by the cooling rate, and the fastest one defined on the mean 
free time scale. This last time scale is the central difference between the two hydrodynamic 
descriptions being discussed here. 




(25) 




(26) 
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C. Comparison of hydrodynamic descriptions 



To compare and contrast the above two hydrodynamic descriptions, it is sufficient to 
consider the homogeneous hmit, for which the details of the constitutive equations are not 
required. Consider ffist the case of zero flow velocities, in which case Eq. (fT5]l gives 

nT = n{r^ + n2T2. (27) 

Equations (fT9|) -( pT]) reduce to constant species densities, and a time dependent temperature 
obtained from 

dtT = -C{T)T, (28) 

while the two temperature description, Eqs. f l22|) - fl24|) . gives also constant species density 
and 

dtTi = -Ci (Ti, T2) Ti, dtT2 = -C2 (Ti, T2) T2. (29) 

Equation fl25]) follows from Eqs. fl29]) via the exact relationships (1251) and fl27|) for the uniform 
case. In this respect, the two temperature description is more complete. However, it is also 
more complex and masks the simplicity of the long time dynamics of ( |28|) by superposing 
on it a higher frequency dynamics, as described in the following. 

An approximate representation for the cooling rates is given in Appendix [B] in the form 

Q 

O = ^(T2)e^(0), (30) 

where 

is a measure of the temperature ratio, and z/ > is some average collision frequency. The 
explicit expressions for u, ^1, and ^2 are given in Appendix El Equations ( l29l) lead to an 
equation for 0, 

9iln0 = -z/(r2)(ei(0)-6(0)). (32) 
There is a stationary solution (f)h determined by 

ei (0/.) = 6 (00 • (33) 

It follows from Eqs. and dHS]) in Appendix M that ^1 (0) - 6 (0) > for > 0,, 

and (0) — ^2 (0) < for < so the stationary solution is approached in general for 
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times greater than the mean free time. This solution is the homogeneous coohng state for a 
granular mixture. 

Stationarity of implies that the cooling rates of each species are the same, 

dtTi dtT2 

^ = ^ ^ Ci = C2 = C- 

J-1 J-2 

The equivalence to C, follows from the relationship nT = riiTi + n2T2. Thus all three 
temperatures are different in general but their cooling rates are the same. The equivalence 
of the cooling rates means only one temperature field is needed to describe the dynamics on 
this time scale, for arbitrary initial preparation. For example, in the two temperature case 
considered above, the homogeneous solution corresponding to (ITS!) is 

- (Ci - C2) <Pw- - (T— + l]p, = 0, (34) 



dT 

where a change of variables has been made from {Ti,T2} to {T, </>}. The solution of this 
equation depends on two time dependent fields, ph = ph{T,T(t), (pit)). But on the mean 
free time scale -> (ph, Ci {4>h) = C2 {4>h) , and so ph{T,T{t),(p{t)) -> ph{T,T(t),(ph). The 
latter is the homogeneous solution to ( [T5|) . the normal solution for the single temperature 
case. 

This analysis of the homogeneous state illustrates the behavior for more general inho- 
mogeneous states as well. The slow dynamics on the scale of the uniformity parameter in 
f fT9|) -f l2T]) is present as well in the two temperature formulation f l22|) - fl24|) . but the latter su- 
perposes on that either initial transients or a continual modulation at the collision frequency. 
Since it is only the slow dynamics that is relevant experimentally, and for justification of 
the normal solution, the two temperature formulation is unnecessarily complex. This does 
not mean that the slow components of the species temperatures are irrelevant. These are 
provided by the single temperature formulation and the normal solution that gives these 
species temperatures as functions of the slow global temperature, Ti(t) = Ti(T(t)). This 



picture has been confirmed by molecular dynamics simulations in the homogeneous case [15 1 
where the rapid approach of T1/T2 to a constant different from unity is observed on the scale 
of the mean free time. Similar results have been seen for an inhomogeneous vibrated mixture 



of inelastic hard disks 16|, showing that the two species temperatures in the steady state 
are determined by the densities and global temperature profiles, being independent from the 
details of the heating mechanism, aside from a boundary layer next to the vibrating wall. 
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Finally, consider the species velocities, described by (123|) . It is clear that the source 
in these equations drives the species velocities towards the common flow fleld u, on the 
mean free time scale. This is similar to the equilibration of the species temperatures, and is 
the same for both molecular and granular systems. As with the species temperatures, this 
extended description yields transients and modulation of the slow dynamics. Instead the 
single flow hydrodynamics describes only the relevant slow components of the species flow 
flelds, through the normal solution, in the form Uj = Uj(T(t), u(t)). 

In summary, the two temperature, two velocity hydrodynamic description of a binary 
mixture subsumes the one temperature, one velocity formulation by providing a more de- 
tailed description of the fluid. However, the additional details are on the short time scale of 
the mean free time where the dynamics is entangled with many other microscopic degrees of 
freedom, and for which the underlying normal solution is not expected to apply. Of course, 
it has the additional computational difficulty of requiring the determination of many more 
transport coefficients in the constitutive equations [l7| , and solving for ten coupled flelds in 
contrast to six flelds in the simple hydrodynamic description. 



V. DISCUSSION 



The hydrodynamic flelds should include all those flelds with the longest time scale, flxed 
by the uniformity parameter. In addition, flelds with shorter times scales may need to be 
included. For example, systems with long internal relaxation times (viscoelastic liquids) can 
be described by non-local in time constitutive equations or, equivalently, by local equations 
with additional flelds. For some non-equilibrium states with inherent dynamics persisting 
in the homogeneous state, such as for granular fluids, non-conserved flelds may be essential. 
Non-conserved species flelds may be useful for cases where their equilibration times in the 
homogeneous state are long, such as for electron-ion two temperature plasmas. 

A more difficult question is whether to include non-conserved flelds whose time scales 
include the mean free time, since many other non-hydrodynamic degrees of freedom are also 
active on this time scale. There does not seem to be any formal objection to including such 
flelds, since the Chapman- Enskog construction of a corresponding normal solution, outlined 
in Sec. HIl can be implemented. However, the normal solution in such cases is still restricted 
to apply only for time scales large compared to the mean free time. Consequently, the 
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new short time scale component in the hydrodynamics is not relevant and can obscure the 
relevant slow dynamics inherent in all fields. 

To make this point in a different way, let {Va} be the phase functions whose average 
values in some nonequilibrium state are the fields {i/a}- It is possible to write the equations 
of motion for {Va} as formally exact Langevin equations {2, 3|, generalizing Eq. ([T]) 

dtVa (r, t) - M„(r, t\ {Vf,}) = (r, t) . (35) 

The functional is similar to in ([T]), except the former is a time dependent functional, 
depending on time explicitly as well as through its functional argument. The sources fa (r, t) 
represent the dynamics of all other degrees of freedom, beyond the set {Va} ■ The solution 
to fl35l) is equivalent to that from Newton's equations, but the dynamics now has been 
separated into two parts that are supposed to be distinct: a slow dynamics associated with 
the homogeneous equation, modulated by a fast dynamics due to fa (r,^)- The average of 
these sources vanish, so the average of Eq. resembles a hydrodynamic description, 

dtVa (r, t) - M„(r, t I {y(s}) . . . = 0. (36) 

The dots denote additional terms due to fiuctuations in the fields, for instance proportional 
to (^(jja — ya)'^), which occur since Ma is a nonlinear functional. Although the sources no 
longer appear their effects are still present in the explicit time dependence of Ma- Only 
after taking time scales long compared to that of the sources does this time dependence 
become negligible, and a hydrodynamic description obtained. This shows that general sets 
of fields can be chosen, but their hydrodynamics is restricted to long times when a simpler 
description may be available in terms of fewer fields. 
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Appendix A: Momentum and energy sources 



The sources in the momentum and energy balance equations are determined in this 
Appendix, and properties (|25|) and ( 126|) are verified. It is assumed that the particles are 
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hard spheres, although the analysis applies as well for general pairwise additive forces. The 
Liouville operator of the system has the form 

i=l i=l j^i 1=1 jj^i i=l i=l j=l 

where qn {Qit) and pn {Pi2) are the position and momentum of particle i of species 1 (2), 
and A^^i and N2 are the total number of particles of each species. The operators Ti^ji = Tju^ 
describe the binary collision between particle i of species k and particle j of species /, 
k,l = 1,2. 



1. Momentum balance equation 

The average momentum density for species 1 is 

Pi(r) = (^5^Pa5(r-q«i)^, (A2) 
where the brackets denote ensemble average. The time derivative is 

dtY>i (r) = -m^^V ■ ( ^ PiiPiiS (r - qa) ) + ( X] 5Z '^'^^^^ ~ 

\ «=1 / \ j^i i=l / 

/ Ni N2 \ 

+ \Y.Y1 ^^^^^^ - ^^1) ) ■ (^3) 
\j=i j=i / 

The first term on the right hand side of Eq. ( 1A3I) gives the gradient of the kinetic part of the 
momentum flux for species 1, while the second term describes collisions among particles of 
species 1. Finally, the last term represents collisions between pairs of particles of different 
species. The second term can be written as the gradient of a flux, 

XI 5Z ^'^^^^ ~ / = 2 ( ^ ^ ~ ^iiiiPii + ^ (r - qji) Tjuipji) \ 

\ j-ti i=i I \ jj^i 1=1 I 

i=l I 

In the last line use has been made of TnjiPji = —Tnjipn (Newton's third law). Finally, the 
difference of the delta functions can be written as a divergence using the identity 

d 

6{r- qn) -6{r - qji) = J dx — 6{r- xqn - (1 - x)qji) 

= -V ■ (q,;i - q^i) / dxS{r- xqn - (1 - x)qji) . (A5) 
Jo 
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Use of Eq. (1A5|) into Eq. (lA4p identifies tlie second term of (1A3P as minus the gradient of 
the momentum flux tensor for species 1. 

A similar analysis of the third term does not lead to a gradient because exchanging i,j 
in the sums leads to two different pairs, for which Newton's third law does not apply. Hence 
this last term is the momentum source in Eq. (1231) . 

/ Ni \ JV2 

Ai (r) = - / ^ 5 (r - q,i) Fa-2 ) , Fa,2 = J2 Tiij2Pii. (A6) 
\ i=i I j=i 

The property follows directly from this result since 

/ Ni N2 

Ai (r) + A2 (r) = - ( ^ ^ (r - qii) Tiij2Pii + S {r - q^^) Tj2iiPj2) 

\ i=l j=l 

/ Ni N2 \ 

= - ( E E - - ^ - ^^'•2)) ^*y2P.i ) . (A7) 
\j=l j=l I 

where momentum conservation, Tjij2Pj2 = —Tiij2Pii, has been used. Using again the iden- 
tity (lASp this becomes the gradient of the momentum flux tensor due to collisions between 
particles of different species. This confirms Eq. fl2^ . 

2. Energy balance equation 

The average energy density for species 1 is 

/ ^1 2 \ 

^iw=(E|;^^(r-q^i)) ■ (A8) 

Its time derivative is given by 

\i=i ^ / i=i ^ 

\ i=i j=i 1 / 

The first term on the right is the divergence of the kinetic part of the energy flux for species 
1. The second term is due to collisions among species 1 particles while the last term is due 
to collisions between species 1 and 2. The second term can be analyzed as above for the 
momentum density, 

A^i A^i 2 \ 1 / 2 7)2 \ 

E E c - *o ) = 2 (E E (m- 1") + Mr - T,,^) 
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\ jjLi i=l I \ j^i i=l 



,2 



(5 (r - Qji) - 5 (r - q^)) 



(AlO) 

In the second equality, use has been made again of Tnji = Tjm and also of 

where Anji is the energy loss by the pair of particles i,j in the collision. The difference 
between the delta functions in the last term of Eq. (lAlOp can be transformed to the diver- 
gence of the energy flux using again Eq. (lASp . Therefore, the second term on the right side 
of (IA9p is the sum of the energy loss due to inelasticity and the energy flux for particles of 
species 1 alone. 

The last term of flAQp cannot be represented as the divergence of a flux since interchanging 
i and j in the summation introduces a different pair. The total energy source for species 1 
is then identified as 

i=i / \i=i I j=i 1 

(A12) 

The cooling rates in the temperature equation are proportional to Wj, (r) 

CiTi = -^w^. (A13) 



The sum of the two sources is 

Afi Ni N2 N2 

Wi 



-j^ / Afi Ni N2 N2 ^ 

:r) + W2ir) = ( E E ^ - ^^1^1 + E E ^ - ^^2) ^^2,2 



2 , 

jj^l 1=1 j^l 1=1 
Using ( lAllI) this becomes 

/ ^ Ni Ni ^ N2 N2 

Wi{t)+W2{y) = - ( 2 E E ~ ^Jl) + 2 5^ ~ ^^2) Ai2j2 

\ j^j i=l jyti 1=1 

Ni N2 

+ E E ^ ~ qj-2) Aiijs 

/ N2 2 

+ (^E E - q^i) - ^ (r - q.-2)) Ta,2^^ ) . (A15) 
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The first term is the total energy loss due to all inelastic collisions, while the last term 
becomes (using flASP ) the contribution to the divergence of the energy flux from collisions 
between particles of different species. This confirms Eq. (^5^. 



Appendix B: Energy exchange rates 

The "cooling" rates for a binary mixture, Q, appearing in (129|) have been estimated from 
an approximate two particle reduced distribution function (local equilibrium or information 
entropy distribution), to give their dependence on the fields nj,Tj and the mechanical prop- 
erties of the particles. It is assumed that the particles are hard, smooth, inelastic spheres 
with masses, diameters, and densities mj,(jj, and rij respectively. The inelasticity is mea- 
sured by restitution coefficients aij for collisions between particles of type i and j, with 
< aij < 1 and a^j = 1 representing elastic collisions. The cooling rates are proportional 
to an average collision frequency u with the approximate results [li 



47rm2 2 / 8^2 



Ci = Ki, X2ina2i\ (l + a2i), (Bl) 

6m V 7rm2 

6 = (1 - «?i)xi/3i 70 + X2^^^ f 1 - /i + 00 - (B2) 

1 + 00 V 4>J 

6 = (1 - «L) ^2/32 + x,^^^ (i + ^ + 0V1 - ^^)■ m 

1 + 00 V ^^ J 

The variable is a measure of the temperature ratio defined in Eq. fl3T|) . a2i = (o"i + o"2)/2, 
mequivnii + m2, Xij is the pair correlation function at contact of a particle of species i and 
a particle of species j, Xi = rii/n, and the constants (3i, (32, 0o, and fi are given by 

n l+0OXll/f^lV 1 + 00X22/^X2^ 

4v2/i X21 Vo'21/ 4v2/i X21 V'^21/ 

1 - a2i m2 , . 

00 = — , n= — . (B5) 

1 + 021 m 

As noted in the text, Q has two contributions, one that is positive and vanishes as (1 — afj) 
in the elastic limit, and another that remains finite in that limit. 
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